Magnetism in SQUIDs at Millikelvin Temperatures 
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We have characterized the temperature dependence of the flux threading dc SQUIDs cooled to 
millikelvin temperatures. The flux increases as 1/T as temperature is lowered; moreover, the flux 
change is proportional to the density of trapped vortices. The data is compatible with the thermal 
polarization of surface spins in the trapped fields of the vortices. In the absence of trapped flux, we 
observe evidence of spin-glass freezing at low temperature. These results suggest an explanation for 
the "universal" 1// flux noise in SQUIDs and superconducting qubits. 
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Superconducting integrated circuits are a leading can- 
didate for scalable quantum information processing [l[. 
A key qubit figure of merit is the dephasing time, the 
time over which a superposition of the qubit |0) and |1) 
states maintains phase coherence. Dephasing is governed 
by low-frequency noise in the bias parameters that con- 
trol the energy separation between the qubit states 0, El • 
In the Josephson flux qubit, the gaussian decay envelope 
of qubit Ramsey fringes is compatible with a magnetic 
flux noise with 1// spectrum and a magnitude at 1 Hz 
of 1 /if&o/Hz 1 / 2 , where <f>o — h/2e is the magnetic flux 
quantum @, [B[. In the Josephson phase qubit, a recent 
experiment used the resonant response of the qubit to 
measure the power spectral density of magnetic flux noise 
directly; again the spectrum was 1 //, with a noise magni- 
tude comparable to that seen in the flux qubit . While 
these experiments identify flux noise as a dominant de- 
phasing mechanism, they offer no clue to its microscopic 
origin. 

The flux noise inferred from recent qubit experiments 
is consistent with the flux noise observed more than 20 
years ago in a series of measurements performed on dc Su- 
perconducting QUantum Interference Devices (SQUIDs) 
cooled to millikelvin temperatures . The noise was ob- 
served to be "universal" , that is, only weakly dependent 
on a wide range of parameters such as superconducting 
materials, SQUID loop geometry and inductance, and 
temperature. While these experiments ruled out many 
potential sources of noise, the microscopic mechanism 
was never identified. There has been recent theoretical 
interest in the possibility that the 1// flux noise is due 
to unpaired spins on the surfaces of the superconductors 
d, 0. Models for 1/f flux noise from surface spins are 
attractive, as they yield a noise power that is only weakly 
dependent on the overall scale of the device, compatible 
with the "universal" character of the noise [6(. To date, 
however, there has been no direct experimental evidence 
for magnetism in superconducting circuits. 

Here we present clear evidence for a high density of 



unpaired surface spins in thin-film SQUIDs. While fluc- 
tuations of these spins produce a small but measurable 
effect, the coherent magnetization of the spins couples a 
very large flux of order 1 $o to the SQUID. Moreover, 
simple measurements of the temperature dependence of 
the magnetization allow determination of the surface den- 
sity of unpaired spins, and indicate that interactions be- 
tween spins are significant. 

We first characterized the flux noise of a two- 
stage SQUID amplifier at millikelvin temperatures. A 
schematic of the experiment is shown in Fig. 1. Both 
SQUIDs are made in the conventional Ketchen-Jaycox 
square- washer geometry (see inset) [ltj. The 110 pH Al 
SQUID is biased with a voltage, and the current through 
the SQUID is monitored with the 350 pH Nb SQUID, 



which is operated in a flux- locked loop ll| . The current 
through the first-stage Al SQUID oscillates with applied 
flux, with periodicity <I>o. When the first-stage SQUID 
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FIG. 1: (a) Schematic of two-stage SQUID amplifier. The 
first-stage Al SQUID is biased with a voltage, and the current 
through this device is read out with a Nb SQUID operated 
in a flux- locked loop. Both devices are operated inside a su- 
perconducting aluminum shield (dashed box), (b) Layout of 
the thin-film SQUIDs. The Al SQUID has inner and outer di- 
mensions 60 and 300 fim; the Nb SQUID has inner and outer 
dimensions 200 fj,m and 1 mm. (c) Power spectral density of 
flux noise of the two-SQUID amplifier measured at 100 mK. 
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FIG. 2: Temperature dependence of flux threading the 
SQUIDs of the amplifier described in Fig. 1. Inset shows the 
I — 5> characteristic of the Al first-stage SQUID measured at 
three temperatures; flux "I>i through the Al first-stage SQUID 
is plotted on the horizontal axis, and current / through the 
Al first-stage SQUID is expressed as a flux $2 coupled to the 
Nb second-stage SQUID and plotted on the vertical axis. 



is biased at a flux (n ± l/4)<I>o, where n is an integer, 
sensitivity to applied flux is maximum. 

The power spectral density of flux noise in the Al 
first-stage SQUID is shown in Fig. 1. The noise scales 
as 1// at low frequencies, with a noise magnitude of 3 
/i'&o/Hz 1 / 2 at 1 Hz. The flux noise is roughly indepen- 
dent of temperature over the range from 30 mK to 500 
mK (not shown) . These results are compatible with pre- 
vious measurements of SQUIDs and qubits @, [H, [fl Q • 

At the lowest temperatures, however, we observe a sur- 
prising dependence of SQUID flux on bath temperature. 
We characterize the temperature dependence of the flux 
in the SQUIDs by tracing out the full current-flux (7— $) 
characteristic of the first-stage SQUID as bath tempera- 
ture is varied. By following shifts in the position of the 
extrema of the 7 — $ curve along both the $ and 7 axes, 
we monitor changes in the quasistatic fluxes threading 
the first- and second-stage SQUIDs, respectively. The 
results of this experiment are shown in Fig. 2. 

The fluxes threading the SQUIDs change dramatically 
as bath temperature is lowered. At the lowest accessible 
temperatures ~ 90 mK, the temperature-to-flux transfer 
coefficients are in excess of 10 $o/K. This behavior has 
been observed in both Al and Nb devices, made both with 
and without a wiring dielectric, and prepared in different 
facilities according to different fabrication recipes. 

Temperature-induced flux drift in SQUIDs has been 
studied in the past. There have been investigations of 
flux shifts due to thermally-driven motion of magnetic 
flux 
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13l |. and due to temperature-dependent critical 
currents in asymmetric SQUIDs [12]. However, we do not 
expect these mechanisms to play a role here, since our ex- 
periments are conducted far below the superconducting 
transition temperature T c . Indeed, the large magnitude 



and low energy scale of the effect observed here suggest a 
new mechanism. Over a broad temperature range below 
500 mK, the flux threading the SQUIDs scales inversely 
with temperature. Paramagnetic impurities in the ma- 
terials of the SQUID would naturally give rise to such a 
signature, and we interpret the 1/T dependence of the 
flux as strong evidence for unpaired spins, most likely in 
the native oxides of the superconductors. However, as 
these experiments are performed in a nominal zero field, 
the source of the orienting magnetic field is not immedi- 
ately clear. 

In order to clarify the source of the temperature- 
dependent flux, we have performed a series of field-cool 
experiments in which a magnetic field Bf c is applied to a 
350 pH square-washer Nb SQUID as it is cooled through 
T c ; the field freezes magnetic flux vortices into the Nb 
film, with density a v ss 7?/ c /<j> [14| . When the device 
is well below T c , the field is removed, and the SQUID 
is maintained in a flux-locked loop as it is cooled to 
millikelvin temperatures. In Fig. 3a we plot the flux 
through the SQUID as a function of temperature for 
eight different values of the cooling field. The cooling 
field strongly affects the temperature-dependent flux, en- 
hancing or even reversing the polarity of the observed 
signal. In Fig. 3b we plot the flux change on cooling 
from 500 mK to 100 mK as a function of the cooling 
field; a linear fit to the data yields a slope of 1.3 $ / m T. 
Clearly, vortices contribute significantly to the measured 
temperature-induced flux shift. 

The linear dependence of flux drift on vortex density 
suggests the following interpretation of the data. As the 
superconducting films are cooled through T c in an applied 
field, vortices nucleate and relatively large magnetic fields 
of order 10 mT are frozen into the films. The observed 
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FIG. 3: (a) Temperature dependence of the flux threading 
a 350 pH Nb/AIOx/Nb SQUID, for different fields B fc ap- 
plied as the device was cooled through the superconduct- 
ing transition (from bottom to top, Bf c — -500, -20, 0, 20, 
50, 100, 200, 500 fiT). (b) Temperature-induced flux change 
A<E> = $(100mK) - $(500mK) on cooling from 500 mK to 
100 mK vs. cooling field Bf c . A linear fit to the data yields 
a slope A$/B /c = 1.3 <E> /mT. 
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signal is due to the magnetization of unpaired electron 
spins on the surface of the superconductor in the strong 
fields produced by the trapped vortices. Indeed, careful 
analysis of the field-cool data allows us to extract the 
surface density a s of spins, a key parameter in models 
of l/f noise from surface magnetism [la ]. We note that 
there is negligible direct coupling to the SQUID from 
spins polarized out of the plane of the superconducting 
films; this is easily understood from reciprocity, as surface 
magnetic fields due to currents in the SQUID have van- 
ishing perpendicular component. However, the topology 
of our device is quite different from that of a continu- 
ous superconducting washer, due to the presence of the 
vortices. The thermal polarization of unpaired surface 
spins in the vortex forces a redistribution of the circulat- 
ing supercurrents in the vortex, due to the requirement 
to conserve magnetic flux. These currents, in turn, cou- 
ple strongly to the SQUID loop. Calculation of the spin 
density from the data of Fig. 3 therefore proceeds in two 
stages: (1) calculation of the coupling between a vortex 
and the SQUID loop, and (2) calculation of the coupling 
between surface spins and a vortex. (A detailed discus- 
sion of this calculation can be found in [l5|). 

A single vortex trapped in the SQUID washer induces 
an effective flux offset in the SQUID; the magnitude of 
the flux offset approaches <i>o and 1 <f>o f° r vortices 
positioned at the outer and inner edges of the SQUID 
washer, respectively. Numerical calculations for a circu- 
lar washer show that this effective flux offset averaged 
over the area of the washer is 0.5 $o for a thin washer 
(with width much less than radius), as expected; for our 
device dimensions (inner radius 100 /jm and outer radius 
500 /xm), we find an average flux offset of 0.14 $o- The 
physical source of the flux is vortex currents that drop off 
slowly with distance from the core and circulate around 
the hole in the SQUID washer. Polarization of the sur- 
face spins couples a flux $„(T) to the vortex, requiring 
the circulating currents to change due to the quantiza- 
tion condition. The change in flux coupled to the SQUID 
is thus $ = -0.14$„(T). 

To calculate $„(T), we rely on reciprocity: knowledge 
of the current distribution in the vortex allows us to de- 
termine the coupling to a spin at an arbitrary location. 
We model the vortex as a hole in a superconducting 
groundplane with radius equal to the estimated coher- 
ence length £ = 30 nm and a thin-film penetration depth 
A = 100 nm. Numerical solution of the London equa- 
tions allows us to solve for the currents in the vortex, 
from which we determine the in-plane and out-of-plane 
magnetic fields B r {r) and B z (r). One can show that the 
flux coupled to the vortex by the spins is given by 



have defined an effective spin polarization P e ff(T): 



$ v (T)=iiBcr s L v P eff (T), 



(1) 
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where B(r) = [B r (r) 2 + -B z (r) 2 ] 1 2 . A change in tem- 
perature yields a change in effective spin polarization 
P e f f (T) in the vortex, which in turn induces a change in 
flux & v (T) coupled to each vortex. The total flux change 
A$ at the SQUID is obtained by summing over all vor- 
tices, taking into account the coupling factor 0.14. Thus, 
we can relate the measured temperature-dependent flux 
to spin density as follows: 



= . U -^i iB a s L v AP eff , 



(3) 



where L„ is the vortex self-inductance, and where we 



where Asq is the area of the SQUID washer. For the 
materials parameters given above, we find AP e ff = 
P ef f (100 mK) - P eff (500 mK) = 0.037 and L v = 0.24 
pH. Using A S q = 0.96 mm 2 and the slope A$/_Bj c = 
1.3 $o/mT from Fig. 3b, we find a spin density a s — 
5.0 x 10 17 m~ 2 . This density of surface spins is compati- 
ble with densities considered in recent theoretical models 
of 1// flux noise from surface spins [1, Q, and with the 
surface spin density postulated in a recent model for spin- 
flip scattering that has been proposed [l(| to account for 
the magnetic-field enhancement of the critical current of 
superconducting nanowires [l7j ]. 

Yet while spin densities of this order of magnitude are 
clearly required to explain the measured l/f flux noise 
and qubit dephasing times in terms of surface magnetism, 
the microscopic physics that drives the spin fluctuations 
is by no means obvious. A recent theory explains flux 
noise in terms of spin diffusion and interaction between 
surface spins via the RKKY mechanism [ll|. A key 
parameter in this model is the characteristic energy of 
spin interactions. In the field-cool experiments described 
above, any possible spin interactions are obscured by 
strong coupling of the spins to the fields of the vortices. 
To probe interactions between surface spins, it is neces- 
sary to examine devices that have no trapped flux. To 
this end, we have performed additional experiments us- 
ing SQUIDs made from narrow (2 /xm) superconducting 
lines. It is energetically unfavorable for a vortex to nu- 
cleate in a superconducting strip of width w provided 
the strip is cooled through T c in a magnetic field that is 
less than the threshold field B thr ~ <S>o/w 2 0, El, HI . 
As the threshold field for our narrow-linewidth SQUIDs 
is around 500 /iT, we expect no trapped vortices when 
these devices arc cooled in a nominal zero field. 

In Fig. 4 we show typical flux vs. temperature data 
for two narrow linewidth Nb/AlOx/Nb SQUIDs con- 
figured as square loops with dimensions d = 50 /xm 
and 200 /im (see inset). The narrow devices show a 
strongly temperature-dependent flux even in the absence 
of trapped vortices. Although the exact magnitude of 
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in the flux vs. temperature data of the 870 pH SQUID 
This is reminiscent of the cusp in the first-order suscep- 
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FIG. 4: Temperature dependence of the flux threading narrow 
(2 fim) linewidth Nb/AlOx/Nb SQUIDs with inductances 160 
pH and 870 pH. Device geometry is shown in the inset. 



the total flux change is seen to vary from cooldown to 
cooldown, suggesting a dependence of the magnetization 
on the details of the cooldown, we consistently observe 
a flux shift of order 10 m$o on cooling to 30 mK. We 
take this as confirmation that the flux drift observed in 
our experiments is not due to thermally active vortices. 
Moreover, these measurements allow us to extract an en- 
ergy scale of spin interactions, as we explain below. 

First, the magnitude of the temperature-dependent 
flux in these devices indicates a high degree of spin or- 
der. Following the analysis of [&], one can show that the 
maximum flux coupled to the SQUID from surface spins 
(achieved when the spins are fully polarized and oriented 
radially) is given by 
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(4) 



The flux scales linearly with device dimension, and is in- 
dependent of the width of the superconducting lines for 
w « d. For the 160 pH SQUID with d = 50 fim, a 
spin density a s =5x 10 17 mT 2 yields $ raal = 280 m$ ; 
the measured flux change of 19 m<I>o therefore suggests a 
high degree of spin polarization at 30 mK. While surface 
spin densities in the narrow SQUIDs could be higher than 
that inferred from the field-cool experiments (the narrow 
devices are not passivated by a wiring dielectric), the 
large temperature-induced flux shift in these devices in 
the absence of a strong applied field suggests an energy 
scale for spin interactions of order 10 mK. Pure dipo- 
lar coupling at the spin densities considered here would 
yield a much lower energy scale, of order 100 fiK. On 
the other hand, the RKKY model of predicts a spin 
interaction strength of 20 mK for the spin densities and 
materials parameters considered here, assuming a uni- 
form distribution of spins through the native oxide; for 
spins concentrated at the metal-insulator interface, the 
spin interaction strength approaches 500 mK. 

Moreover, we note the appearance of a cusp at 55 mK 



tibility of a spin glass at the freezing temperature [21 1. 
This feature again indicates the importance of spin in- 
teractions at millikelvin temperatures and suggests the 
possibility of spin-glass freezing of the unpaired surface 
spins on the timescales of our experiment. 

In summary, we have demonstrated that superconduct- 
ing circuits are magnetically active at millikelvin temper- 
atures. This observation points to a microscopic explana- 
tion for the excess 1/ / flux noise in Josephson circuits, a 
dominant source of dephasing in superconducting qubits 
and an open question for more than 20 years. We observe 
evidence of interactions between surface spins and pos- 
sible spin-glass freezing at millikelvin temperatures. A 
high density of surface spins could play a deleterious role 
in other schemes for quantum information processing in 
the solid-state, notably those based on electron spin. It 
is possible that suitable surface treatments of the super- 
conducting films will lower the density of magnetic states, 
leading to superconducting devices with lower noise and 
solid-state qubits with improved coherence times. 
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